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Abstract 

We derive and solve the Hamiltonian flow equations for a Dirac 
particle in an external static potential. The method shows a gen- 
eral procedure for the set up of continuous unitary transformations to 
reduce the Hamiltonian to a quasidiagonal form. 
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The well known Foldy-Wouthuysen (FW) transformation is a unitary 
transformation that separates big and small components in the Dirac equa- 
tion for a particle interacting with an external electromagnetic field. If the 
external field is stationary then the FW transformation reduces to a unitary 
transformation of the Hamiltonian H of the Dirac particle 

H' = e'^{a{p- eA{x)) + (3m + eAo{x))e-'^ (1) 

where 5* can be expanded in a sequence of hermitean operators in powers of 
1/m. Recently a novel method to diagonalize Hamiltonians via continuous 
unitary transformations has been proposed 0, where the Hamiltonian H 
and the S are considered to be functions of a continuous variable, the flow 
parameter /. This new method has been applied to several topics like the 
n— orbital model superconductivity and the electron phonon interaction 
[0, We want to apply it here to the Dirac particle in an external fleld in 
order to explore its possibilities with a well known problem. Thereby we can 
gain experience in solving problems like the elimination of negative energy 
states in near light cone QCD where the solution has yet to be discovered. 
For inflnitesimal steps exp{ri{l)dl) a continuous unitary transformation of the 
Hamiltonian H leads to the so-called Hamiltonian flow equation 

^ = W),H{1)] (2) 

with initial condition H{0) = H . The transformation matrix ri{l) is antiher- 
mitian and related to S as follows 

The parameter / has dimension of energy''^ and characterizes the amount of 
diagonalization which has been achieved. In the energy representation the 
matrix elements of the Hamiltonian H will be grouped in a band of width 
AE = — ^ around the diagonal, if ri{l) is chosen in the correct way. The 

result of the transformation depends critically on the choice of ri{l) as a 
function of H{1). In ref. the following expression has been proposed: 

^{l) = [Hn{l),H{l)] (3) 
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where Hnil) is the diagonal part of H{1). Indeed such a choice of rj{l) leads 
to a quasidiagonalization of H{1) when / — >• oo. From eqs. (0) and (0) we 
obtain that 



Jl 



d 
Jl 



tiHUil) = 2ti{[Hil),HDmHD{l),H{l)]) > 



Thus the sum of the absolute squares of the nondiagonal elements decreases, 
whereas the sum of the diagonal matrix elements squared increases with I. 
However, this line of reasoning only holds if one can solve eqs. (@),® exactly. 
For approximate solutions of eqs. (Q) and divergencies can arise|0]. 

In this report we want to present a continuous transformation of the Dirac 
Hamiltonian via flow equations the derivation of which is guided by the 
FW transformation. The basic idea of the FW transformation is a decom- 
position of the Hamiltonian into two noninteracting parts corresponding to 
positive and negative energies. Therefore, the transformed Hamiltonian must 
commute with the /3 matrix and we choose as ansatz ?7(/) in the form 



The Hamiltonian H{1) can be presented as a sum of an even operator 
£{l) and odd operator 0{l), where the even or oddness is defined by the 
commutation relations of the respective operators with the f3— matrix. 



v{l) = [Pm,Hil)] 



(4) 



H{1) = 0{l)+£{l), 
£{1)I3 = +p£{l), 
0{l)f3 = -(30{l). 



The initial conditions for the evolution equation 




are obtained by using the plus commutator [a, /?] 



+ 



0. 



£{0) 
0(0) 



(y.{p — eA). 



Then our ansatz for the transformation matrix becomes 
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The flow equation can be split up into the following system of two equations 

4mf30\l), (5) 
2mP[Oil),£il)]. (6) 



dl 
dO{l) 



dl 

As we will show now, the system of eqs. (|^),(|) can be solved perturbatively 
in 1/m. It is convenient to introduce a dimensionless flow parameter A = 
I * rn?. Since £^(0) = /3m + cAq the expansion of £{\)/m in a series in 1/m 
contains terms starting with the zeroth order term 

-£{\) = £o(A) + -^i(A) + ^^2(A) + . . . , (7) 
m m m^ 

whereas the expansion of 0{\)/m starts with the first order 

-0{\) = -Oi(A) + ^02(A) + . . . (8) 
m m 

Due to eqs. (0) and (||) we have in n-th order : 
d 

—Sn{\) = 4/?5^0fc(A)a-;t(A) (9) 
k=i 

-^a(A) = -AOn[\) + 2(3Y.[Ok{\),£n-k{\)] (10) 

"'^ fc=l 

The solution of these equations are 



£n{\) = £n{0)+A(3[ dA'(^Ofc(A')0„-fc(A')) (11) 

k=i 

a(A) = a(0)e-^^ + 2/5e-^M rfA'(^ [e^^ 0,(A'), ^„-fc(A')]) (12) 



with initial conditions 



^o(O) = /3, £i{0) = eAoix), £„(0) = if n > 2, 
Ci(0) = a(p - A(f)), C„(0) = if n > 2. 
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One sees that On(A) exponentially goes to zero when A oo. By virtue 
of this behavior integrals in eq.(|lT]) do not diverge and £^n(A) is finite when 
A — » oo. 

For the first four orders we have 

^o(oo) = So{0)=(3, 

^i(oo) = fi(0) = eAo(f), 

^2(00) = ^(3Ol{0) = ^(3[{p-eA{x))'-eaB{x)], 

£,{oo) = ^[[Oi(0),^i(0)],Oi(0)] = 

= __a^(v X E{x)) - -a{E{x) X {p- eA{x))) - -iyE{x)). 

These terms reproduce well known Foldy-Wouthuysen Hamiltonian The 
formulae (pl]) , ([l2|) let us simply get expressions in any order in 1/m. 

Similarity to the renormalization group equations the Hamiltonian fiow 
equations approach the correct solution in infinitesimal steps. We have de- 
rived in this new method a solution of the fiow equations (H) and @ for the 
problem of a Dirac particle in an external potential. This solution has a more 
general validity. We have shown how to choose the diagonalizing matrix 77, 
if the Hamiltonian can be represented as a sum of even and odd operators 
under commutation with some operator D. Choosing the operator r]{l) in the 
form [D, H{Vj\ allows us to transform the Hamiltonian if to a form which 
commutes with D whereas the anticommuting part goes to zero when the 
fiow parameter goes to infinity. 
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